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Dynamic Peierls-Nabarro equations for elastically isotropic crystals 
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The dynamic generalization of the Peierls-Nabarro equation for dislocations cores in an isotropic elastic 
medium is derived for screw, and edge dislocations of the "glide" and "climb" type, by means of Mura's eigen- 
strains method. These equations are of the integro-differential type and feature a non-local kernel in space and 
time. The equation for the screw differs by an instantaneous term from a previous attempt by Eshelby. Those 
for both types of edges involve in addition an unusual convolution with the second spatial derivative of the 
displacement jump. As a check, it is shown that these equations correctly reduce, in the stationary limit and 
for all three types of dislocations, to Weertman's equations that extend the static Peierls-Nabarro model to finite 
constant velocities. 

PACS numbers: 61.72.Bb, 61.72.Lk, 62.20.F— 



I. INTRODUCTION 

Plastic deformation in crystals occurs as dislocations move 
through the material under an applied stress. 1,2 Major quan- 
titative progresses in plasticity modeling arose with the out- 
break of the Peierls-Nabarro (PN) integral equation. 3-5 Aimed 
at computing dislocation core shapes, this equation establishes 
a quantitative link between atomic forces, nowadays described 
by means of the material-dependent 7 potential (a "reduced" 
lattice potential specialized to shear deformations) 6 and the 
dislocation core structure. Since, numerous refinements of 
various nature improved the agreement between the PN model 
and molecular statics simulations, though best matches with 
experiment for the core width and the Peierls stress 3 are ob- 
tained so far not by using the PN model, but by addressing 
ab initio the full three-dimensional structure of dislocations 
cores. Still, in spite of known drawbacks, the PN equation 
remains widely used. These questions have recently been re- 
viewed by Schoecki 7 - 

Yet, the dynamic instance of the Peierls-Nabarro equation 
appears as a long-standing elusive issue in dislocation theory. 
To date, simulations (using molecular dynamics^iS or phase- 
field methods 11 ) constitute the privileged path to specific dy- 
namic core-related phenomena. Among the latter are the long- 
hypothesized transonic or supersonic transitions^— first re- 
ported in atomistic simulations 17-20 but only recently exhib- 
ited experimentally (in a two-dimensional plasma crystal). 21 
However, their current understanding is far from complete, 
given the variety of situations that can arise depending on 
the nature of the crystal and environmental conditions. An- 
alytical progresses have also been made in this direction, but 
most often at the price of restrictive approximations for the 
dislocation core. 22 25 It should be noted that the major diffi- 
culty involved in sonic transitions in real crystals, recognized 
long ago, 14 resides in that a dislocation can be part subsonic 
and part supersonic at the same time due to spatial-dispersion 
effects. 19 Tackling this problem is not attempted here. How- 
ever, in the above context and even within a non-dispersive 
approximation, a one-dimensional dynamic PN equation that 
leaves all freedom to the core shape would certainly constitute 



a useful additional tool. 

The truth is, in his classical '53 paper on the dynamic mo- 
tion of dislocations, 13 Eshelby did write down a dynamic gen- 
eralization of the PN equation for the screw dislocation. How- 
ever, acknowledging its complexity he did not use it, focus- 
ing instead on an equation of motion for screw dislocations 
under an assumption of rigid core. A dynamic PN equa- 
tion for edge dislocations was never proposed, and in practice 
the only velocity-dependent PN equations studied so far are 
Weertman's 16 and its modifications, 26 which apply to constant 
velocities only. Despite a number of analytical explorations of 
the dynamic regime, this gap has not been filled in yet. 

Quite unexpectedly, a close examination of Eshelby 's dy- 
namic equation for screws 13 leads one to conclude that it 
does not reduce to Weertman's equation in the stationary 
limit. This can be seen from the calculations in Appendix 
CI. One clue to the reason of this discrepancy is provided 
by the recent observation 27 that classical static expressions 
for dislocation-generated displacements, such as that found 
in Refs.[l]andUmiss one term (a distribution) that represents 
the non-elastically relaxed slip. This term proves irrelevant 
to the standard static PN model, and cannot be spotted from 
the static elastic strains alone, since it preserves registry (Sec. 
Hil l. It is shown below that one term of similar origin is rel- 
evant to dynamic calculations, and provides the explanation 
for the above discrepancy. With this observation, the Green's 
function machinery 2 can safely be harnessed to produce the 
desired dynamic PN equations for screws and edges, that cor- 
rectly admit Weertman's equations as stationary limits, pro- 
vided that attention is paid to distributional parts in carrying 
out various Fourier integrals (Sec.Hllil. 

For convenience, indices i = x,y, and z or 1, 2, and 3 are 
used indifferently hereafter. To ease calculations, most of the 
integrals are read in Ref.H^. Reference is made to these inte- 
grals by their book classification number, preceded by "G.R". 
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n. GREEN'S FUNCTION APPROACH TO DISLOCATIONS 
A. Eigenstrains and dynamic Green's function 

Inclusions or defects such as dislocations produce distor- 
tions in their surrounding medium. The total distortion (3 
is the gradient of the material displacement u, such that 
/3y (x) = djUt. Its symmetric part is the total strain = 
(1/2) (Pij + f3ji). Assuming small deformations, the total 
distortion produced by a defect can be written as the sum of 
a linear elastic distortion /3 e , and of a "nonlinear" part (3* 
usually called eigendistortion^^ @ij = f3fj + f3*j, none of 
the latter quantities being a gradient in general. Whereas the 
eigendistortion represents a purely geometric, rigid, i.e., non- 
elastically relaxed, 30 contribution to the total distortion that 
results from the insertion of the inclusion, the elastic distor- 
tion represents the elastic relaxation correction that confers to 
j3 a gradient character. A similar decomposition holds for the 
strain: e l} = £?•+£■,•. 

The Green's function approach to dislocations 2 consists in 
representing the dislocation by an eigendistortion (localized 
on a the glide plane) whose physical interpretation is given 
below (Sec. MB I) and in computing the induced displacement 
field u using an elementary solution of the equations of elas- 
ticity. The total distortion /3y ensues, from which (3 e is ob- 
tained by subtracting (3* . Finally, the stress a follows from 
the linear-elastic strain e e = sym f3 e and linear elasticity, as 

Uij = Cijkie e kl = Cijki(eu-£*ki) = Cij k i{d k ui~ (5* kl ), (1) 

where Cijki — Cijik — Ckiij are components of the elastic 
tensor. Momentum conservation in the form finery — pdfui, 
where p is the mass density, is written as 

CijkidjdkUi - pd 2 t Ui = djTij, (2) 

where ry = Cijki/3'^. In an infinite medium, the Green's 
function of the displacement, G(x, t), is the solution corre- 
sponding to a point-like source located at the origin of space 
and time (the minus sign is conventional)^ 

C ijk idjd k Gi m (x,t) - pd 2 G m (x 7 t) = -6 im 6(x)5(i). (3) 

The following space-time Fourier transform (FT) conven- 
tion is used hereafter (/ is an arbitrary function): 



/(*,*) 



d A k duj 
(2tt) 3 2^ 



/(k,w)e i(k - x - wt) . 



(4) 



Introducing the acoustic tensor N of components iVjy — 
Cikijkkki and the identity matrix I of components Sij, the so- 
lution to Eq. (0 reads in matrix notation 



G(k,w) = (N — pcu 2 \) 



(5) 



It is convenient for the problem at hand to work in the mixed 
"space Fourier modes/time" representation. By convolution 
of the elementary solution, the solution to Eq. (01 is obtained 



where the integrals run from — oo to +oo, and the stress fol- 
lows from Eq. ([]]). Henceforth, we confine ourselves to the 
simplest isotropic case for which 



Cijki — ^SijSki + p(SikSji + S^Sjk), 



(7) 



where p is the shear modulus and A is the Lame coefficient 
Introduce moreover the shear and longitudinal sound veloci- 
ties cs = 



J pj p and cl = y(X+ 2p)/ p. The inverse in Eq. 
(0 is immediate in the basis of longitudinal and transverse 
projectors with respect to k = k/fc. Thus, 



N = pk 2 (I - kk) + (c L /c s ) 2 kk 



(8) 



and the dynamic Green's function of the displacement reads 



G(k,w) = (N — plo 2 \] 



I kk 



kk 



k 2 - (oj/c s ) 2 c[ k 2 - (w/c L ) 2 



(9) 



Its static limit is more conveniently expressed in terms of the 
Poisson ratio v = A/[2(A + p)], such that c|/c 2 = (1 — 
2u)/[2(l-u)], 



Gy(k) = G ii (k,w = Q) 



1 

pk 2 



2(1- v) 



(10) 



Let 6(x) denote the Heaviside function. Inverting the time- 
Fourier transform in Eq. (O with a suitable choice of con- 
tour in the cj-complex plane 31 yields the following retarded 
Green's function, which describes waves going away from the 
source: 



Gi 3 -(k,i) 



9(t)c 2 



pk 



1 

cs 



sm(cskt)(Si 



ki kj j 



sin(cL kt)kikj 



(11) 



B. Volterra dislocations and importance of history 

The problem of finding the fields associated to a disloca- 
tion with an extended core is most efficiently split up in two 
steps. First, a solution is derived for a Volterra dislocation 
of infinitely narrow core, which only slightly complicates the 
above calculation for a pointlike source. In the second step, a 
convolution product of the obtained elementary solution with 
the shape of the extended core, considered as a superposition 
of Volterra dislocations, is taken according to the superposi- 
tion principle of solutions of linear elasticity. This approach, 
introduced by Eshelby, 5 is well suited to obtaining the Peierls- 
Nabarro integral equation. Indeed in this equation (one of 
stress balance) where the core shape itself is the unknown, the 
convolution integral cannot be explicitly evaluated in general. 

The eigendistortions associated to the three relevant types 
of rectilinear infinite Volterra dislocations, with dislocation 



iii(x,i) 



df 



d 3 fc 



Ji.ik.l /':^-,!k./') f ' k ' 



ikx 



(6) 
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a) Screw b) Glide edge c) Climb edge 



FIG. 1. (color online) Three different types of Volterra dislocations 
(core of null width), as computed using Eqs. il6i , < I23 b and d24t with 
v = 0.3, e a = 1(T 4 and b = 1 (|4 \y\ < 5) (Ref.ll 

)■ 

line along the Oz axis, are represented as follows (the core 
lies at the origin of the Cartesian axes): 

p* j (x.)=b6(y)6(- X )S i3 6 j2 (screw), (12a) 
p* j (x)=bS(y)9(-x)S il 5 j2 (glide edge), (12b) 
P* j (x)=b6(y)9(-x)6 i 25 j 2 (climb edge). (12c) 

The norm of the Burgers vector b is b. For the screw, glide 
edge and climb edge, the nonzero component of the Burgers 
vector is b 3 , b\, and b 2 , respectively. The slip plane, where the 
material displacement u experiences a discontinuity, has been 
chosen as y = in all cases. In Eq. ( fT2l . index i refers to 
the components of the Burgers vector, and index j to the nor- 
mal to the plane of discontinuity. Hence, the Burgers vector of 
the screw is parallel to the dislocation line (in this case, the no- 
tion of slip plane proceeds from usual considerations about ex- 
tended loops 1 ), that of the glide edge is orthogonal to the line 
and contained in the slip plane, whereas that of the climb edge 
is orthogonal to the slip plane. FigureQ]illustrates these three 
types of Volterra dislocations in which the reader familiar with 
fracture theory will recognize dislocation counterparts of the 
three conventional modes of fracture, 33 the correspondence 
being (mode III — tearing, mode II — sliding, and mode I — 
opening) o (screw, glide edge, and climb edge). The two 
kinds of edges are of distinct nature: with our sign and ori- 
entation conventions, the positive glide edge is obtained by 
compressing the half-space y > along Ox in the positive x 
direction, whereas the positive climb edge results from insert- 
ing one extra half-plane of atoms along the negative x semi- 
axis, which requires to "open" the pre-existing lattice in the y 
direction (black region in Fig.QJ). As is shown below, these 
differences express themselves in the analytical expressions 
of the displacement fields. In usual conditions, the climb edge 
moves by "climb", an essentially diffusive mode that involves 
slow migration of atoms and interstitials. In quasi-static situ- 
ations, climb dislocations are constituents of sessile prismatic 
loops or Frank partials. 1 However, Weertman has emphasized 
their potential importance in fast dynamics as well^ This 
question is addressed in Sec. IIIIDI 

Equations ( fT2b are particular cases of the more general ex- 
pression for a dislocation lina 2 - 

/3*(x) = b<g>n(x)5 s (x), (13) 



where 5s is a Dirac distribution localized on the discontinuity 
surface S of normal n(x). 34 The relevance of the location of 
this discontinuity surface to non-uniform motion is now dis- 
cussed. Relative to the pristine crystalline state, the structural 
modification generated by the presence of a dislocation can 
be seen as the cumulative effect of elastic atomic displace- 
ments, produced by the dislocation core and associated to 
long-range stresses, and of permanent (irreversible) displace- 
ments of atoms accompanying dislocation motion from its nu- 
cleation location to its current location. The above eigendis- 
tortions are associated to the permanent displacements. By 
definition, the integral of the relative material ("atomistic") 
displacements, 

i> (14) 

where C is a closed contour surrounding the dislocation, is 
non-zero, equal to the local Burgers vector, and tends to the 
true Burgers vector as the loop radius goes to infinity^ This 
non-zero value materializes a discontinuity of the total ma- 
terial displacement in the medium. From a physical stand- 
point, once the atomic perturbations generated by the dislo- 
cation motion have been damped, crystal integrity is restored 
behind the dislocation, and the permanent displacements are 
not observable. For this reason, the above integral is contour- 
independent for contours large enough compared to the core 
size, and provides no information on the trajectory followed. 
Alternatively, the dislocation can be considered as constitut- 
ing the boundary line of the surface S in Eq. (fT3l . so that 
prescribing this surface removes in practice the indeterminacy 
of the position of the displacement discontinuity. Therefore, 
in the static limit and for an infinitely thin dislocation line, the 
surface can be chosen arbitrarily since the location of the dis- 
location line alone uniquely determines the elastic strains and 
stresses. 

Resolving this arbitrariness by deciding to localize /3* on 
the geometric surface spanned by the dislocation line dur- 
ing its motion puts in information about the trajectory in the 
problem. 2 For instance in expressions (1121 1. the dislocation 
'comes' from x — —oo. While this information is irrelevant 
in the static case for infinitely thin dislocations, this is not true 
any more in fast dynamics, or for dislocations with extended 
core such as those considered in the Peierls-Nabarro model. 
In the first case, the atoms perturbed by a dislocation passing 
by oscillate and act as wave sources as long as the perturba- 
tion is not fully damped; in the second case, the stacking fault 
that constitutes the core locally stores potential energy, in par- 
ticular in case of dissociation and emission of one partial dis- 
location, and acts as a continuous stress source in the zones 
where the eigenstrain varies much. 

C. Static Peierls-Nabarro equation 

For clarity and further reference, the method to obtain the 
static PN equation with the Green's function method is briefly 
reviewed. The total displacement of a screw dislocation, due 
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to Burgers, 35 is given in all reference textbooks (e.g., Refs. 1 
and 12) as 

. . bd> b , . b y 
«3W = — = — arg(x + ly) = — arctan-, (15) 

where <f> is the polar angle in the [x, y) plane. However, with 
the principal determination of the arctangent, the latter expres- 
sion in Cartesian coordinates is incomplete. Imposing a cut on 
the negative x semi-axis, the correct result instead reads^ 

u 3(x) = tt arctan - + - siga(y)6(-x), (16) 
2ir x 2 

where the distributional part represents the permanent dis- 
placement. The relative slip 77 between both sides of the slip 
plane is 

r)(x) = u 3 {x, y = + ) — u 3 {x, y = ()-) = b6(-x). (17) 

Even though Eq. ( fTo*b is pretty obvious form the first equality 
in Eq. ( TT3T > and the above elementary remark, retrieving Eq. 
( fT6] l using the static Green's function proves a useful exercise 
prior to considering edges and dynamic calculations. Indeed, 
the calculation is given by Mura (Ref. 2, p. 17), with again Eq. 
(fT~5b as a result. While the cause here may reside in the tables 
used by this author, the correct calculation is reproduced in 
AppendixlAlfor definiteness. 

The nonzero components of the total distortion are obtained 
by differentiation of Eq. ( IT6l ) : 



b 

2^ 

b 



ttS(x) sign(y) 

x z + y^ 



~ sign(y)<%) 



y 



27r x 2 + y 2 ' 
b x 



2tt x 2 + y 2 



b5(y)9(-x). 



(18a) 
(18b) 



In these expressions, use has been made of the identity 
arctan a; + arctan 1/x — (tt/2) signs, from which follows 
the derivative (arctan 1/x)' = irS(x) — 1/(1 + x 2 ). The 
standard textbook expressions of the elastic strain follow from 
e e = sym(/3jj — f3*j). It should be noted that distributional 
parts cancel out, and are absent from the latter expressions, 
consistently with the fact that a static strain does not depend 
on history (see previous section). Stresses are obtained by 
multiplying the (shear) strain components by 2fi, asi 



[ib y 



fib x 



2ir x 2 + y 2 ' zy 2tt x 2 + y 2 ' 



(19) 



To construct the PN equation, information on the slip plane is 
reintroduced by computing the stress at y — ± . Thus 



<Jzx{x,0 =-— lrm-2- — 2 

2tt 2/^0 x A + y z 2 



Trsign^)^) 



fib x fib 1 

a zv (x, 0) = — lim — = = — — p.v. — . 

VK ' ' 2vr y^a x 2 +y 2 2it v x 



(20a) 
(20b) 



where p.v. stands for the principal value. Given expression 
(TTTb of the differential slip, the stress [Eq. ( l20bH produced by 
the dislocation on its slip plane is rewritten as the convolution 
product 



a zy (x 7 0) = -^-p.v. / dx 
2ir J 



x — x' 



(21) 



This expression now holds for any core shape function rj(x). 
Adding an applied resolved shear stress a a (x) to Eq. (l2Tb . and 
balancing their sum by the 6-periodic pullback force of atomic 
origin which derives from the stacking fault 7 potential, here- 
after denoted by f(rj), the static PN equation for the screw is 
obtained: 



2tt 



p.v. 



Ml¥-L+<j a (x) = f(T,(x)). (22) 

-r^t & X 



At this point, boundary conditions express the history of dislo- 
cation formation; typically rj(— 00) = rfp+b and r?(+oo) = rfo 
for a single dislocation coming from x = —00 (Ref. [261) or 
7^(±oo) = 770 for a dipole, 4,36 where 770 is the homogeneous 
solution such that a a = f'(rjo). 

In the glide edge case, the nonzero displacement compo- 
nents are u x (x,y) and u y (x,y). The Fourier integrals in 
Mura's method are slightly more complicated, but similar to 
that for the screw. One finds 3 -^ 



u x (x,y) 



u y (x,y) 



^- arctan - + ^ sign(y)0(-x), 
2-k x 2 



b 1 xy 

4.7T 1 — v x 2 -\- y 2 

rc 

b l y 2 

An 1 — v x 2 + y 2 



(23a) 



(23b) 



where e g is of order the inverse of half the system size. The 
sole difference between the present approach and classical 
results is the presence of the additional distributional term 
(b/2) sign(y)6(—x) in u x . We recall that because of one 
divergent integral, u y is determined only up to an additive 
constant that blows up as e g —> 0. For this reason, differ- 
ent equivalent forms of u y are found in the literature. 5,3S ' 38 ' 39 
This complication, linked to torsion, is well-documented (see 
Ref.Qlp. 78). Equation d23b| ) is the form obtained by Eshelby 
and Mura. 2,5 Analogous expressions for the climb edge can be 
written as: 



U x {x,y) 



b 1 



4.7T 1 — v x 2 + y 2 



(24a) 



b 1 - 2v 



u y (x,y) 



(24b) 



8tt 1 - 

b 1 xy 

Att 1 — v x 2 + y 2 

b x b . 

- — arctan - + - sign(y), 
2tt y 4 

where e c = e g exp 1/(1 — 2v)£L When rotated clockwise by 
a angle 7r/2, i.e. subjected to substitutions (x, y) — > (—y, x) 
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and (u x ,u y ) —> (—u y ,u x ), Eqs. (l24l i become identical to 
Eqs. d23"l l. up to differences in e c and in the distributional part. 
These differences are mathematical renderings of the differ- 
ent nature of glide and climb edges emphasized in Sec. Ill Bl 
though both lead (in the static case only) to identical elastic 
strains and stresses characteristic of an edge dislocation. 

The ensuing static PN equation for the glide edge, for which 
the resolved stress is a xy , is identical to Eq. ( 1221 save for a 
prefactor 1/(1— v) in front of the integral and for the definition 
of 77(2) = u x (x, + ) — u x (x, CP). The static PN equation for 
the climb edge, driven by a tensile a yy stress component is 
formally the same as that for the glide edge, but with now 
r](x) = u y {x, + ) — u y (x, 0~). However in the latter case the 
relevant specialized lattice potential, linked to the introduction 
of interstitials, has not been properly defined to date to this 
author's knowledge. 

Thus, in all three cases the distributional term plays no part 
because it does not show up in the elastic strain. The situation 
markedly changes in dynamics. 



III. DYNAMIC PEIERLS -N AB ARRO EQUATION 

The dynamic calculation is quite analogous to the above 
procedure, using the dynamic Green's function [Eq. ( TTTT il in- 
stead of the static one. The main difference is that time- 
dependent distributional contributions, which generalize the 
above static ones, are no more irrelevant and provide impor- 
tant contributions to the PN dynamic equation. The difficulty 
mainly resides in computing cumbersome Fourier integrals. 
One approach could consist in using the Cagniard-de Hoop 
method, as proposed by Markenscoff and Clifton 40 to address 
dynamic dislocation problems. However, using "brute force" 
and reference tables of integrals proved more expedient for 
the case at hand. The edge cases require one key integral that 
we could not find in tables, which is computed by means of a 
differential equation. Appendix [B] contains a detailed sketch 
of these calculations. Hereafter, |x|| = x 2 + y 2 . 



A. Principle 

The dynamic PN equation is obtained in the following man- 
ner. First, one computes the elementary stress field produced 
by a time-dependent eigendistortion j3*j (x, t) representing an 
"instantaneous" Volterra dislocation located at x = y = 0, 
and present at t = 0. We take /3j*(x, t) equal to any of Eqs. 
dT2l . multiplied by the Dirac impulse 8(t). Eshelby and more 
recent works 22 instead consider elementary Volterra disloca- 
tions proportional to &(—t), but the present approach simpli- 
fies the calculation of the stresses in the perspective of obtain- 
ing a PN equation. Mutatis mutandis, we then follow Eshelby 5 
by appealing to the identity: 



rj{x,t) = r?(+oo,t) 



(25) 



j dTdx'6(-(x- x'))5(t-T)^(x',T). 



Here and below, double integrals over space and time have im- 
plicit bounds ±00 unless otherwise stated and are denoted by 
a single integral sign. The integral term expresses the spec- 
trum of instantaneous Volterra dislocations associated to the 
core shape function 77. No dislocation is associated to the ho- 
mogeneous slip i](+oo,t). Invoking linear superposition, if 
a elem (x, t) is the shear stress on the slip plane generated by a 
Volterra dislocation b9(—x)5(t), the shear stress generated by 
the continuous slip r](x, t) = Ui(x, y = + ) — Ui(x, y = 0~) 
reads, by Eq. ( 125) , 

a(x,t) = -\ I dTdx'a elem (x-x',t-T)^(x',T). (26) 
J ox 

An applied inhomogeneous stress in the bulk moreover pro- 
duces a stress a a (x,t) on the slip plane, to be added to Eq. 
(1261 1. Balancing the resulting expression by the pull-back 
stress yields the dynamic PN equation, 



-\ [ drdx'a e,em (x - x' ,t - t)^L{x' ,t) 
b I ox 



+ <j a (x,t) = f'(r){x,t)), 



(27) 



where r;(±oo,i) is such that <j a (±oo,t) = f'(r](±oo, t)). 
We now proceed to determine er elem for the different kinds of 
dislocations. 



B. Screw dislocations 

Then, the displacement associated to the instantaneous 
screw reads (see Appendix lB 11 1 



xy 



0(cst-|x| a ) 



(cgf-y 2 ) V^ 2 -^ 



T7sign(y)9{~x)d{c s t - \y\) 



(28) 



In this expression, the Dirac term is the dynamic counterpart 
of the static distributional term in Eq. (fTol i and represents a 
wave leaving the slip plane orthogonally to it. The associated 
elementary shear on the slip plane y = follows as 



a e,em (x,t) = a zv (x,y = 0,t) 



= lim a 
Introducing the kernel 

K(x,t) = 
one directly finds: 



du. 



(x,y,t) - p* Jx,y,t) 



dy 

x 0{c s t-\x\) 
2c s t 2 ^2 _ ~2 ' 



(29) 



(30) 



af™(x,y = 0,t) - ^K(x,t) - ^-6{-x)5' {t). (31) 



7T 

produces 



2c s 



Applying Eq. 

a zy (x,t) - 

+ £- s J dTdx ' (-( x - x ')) S '( t - T )^( x '^)- ( 32 ) 



At dx 1 K(x — x , t — t) — (x', t) 
ox 
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Assuming that dr//dt(+oo,t) — 0, the second integral re- 
duces to 

dr dx' 0(-(x - x'))5\t - t)^(x', t) 
+ °° d 2 n dn 

dx 'idt^ = -^ (33) 

which gives the time-dependent stress 



a zy (x,t) 



dr da;' K(x — x ,t — t) — (x, t) 
ox 



(34) 



Apart from the presence of the driving stress er a , this equa- 
tion differs from Eshelby's [Eq. (21) of Ref. 1 1311 by the in- 
stantaneous term proportional to drj/dt. This term should 
be replaced here and henceforth by drj/dt, where rj(x,t) — 
rj(x,t) — Tj(+oo,t), whenever boundary conditions require 
dr)/dt(+oo,t) ^ 0. 



C. Glide edge dislocation 

Dynamic displacement fields for the instantaneous glide 
edge are derived in Appendix lB 21 They read 



Hence, from Eq. ( |27| ). the dynamic PN equation for the screw 



fi f dij 

— / dr dx K(x — x .t — t) — (x .t) 

■K J OX 

-£- s ^(x,t)+a a (x,t) =/' (r](x,t)). 



(35) 



Ul (x,t) = ^0(t) 

Z7T 



Z7T 



2xy 



^ 2 * 2 -W'.(c Li -|x| 2 ) 



2clt 2 - |x 



xy 



9(c s t-\x\ 2 ) 



6(c s t-\x\ 2 ) 



(4t 2 -y 2 ) v/^ 2 -|x| 



+ wsign(y)d(-x)6(c s t - \y\) 



c s x 2 \x\ 2 -c 2 t 2 (x 2 ~y 2 ) 



CL 



0(c L t - |x| 2 ) 



a; 2 |x|| — c 2 t 2 {x 2 — y 2 



-0(c s i-|x| 2 ) 



(36a) 



(36b) 



The corresponding expressions for the distortions and stresses are easy to compute but lengthy so that only a xy (x, y = 0), the 
relevant stress for the PN equation, is reproduced here. Using a xy — ^{(3% y + (5 yx ) — fJ.{u x . y + u VlX — (3* y ) with (3* y — 
W(-x)6(y)6(t) yields 



al 1 ™(x,y = 0,t) = ^ 

y 7T 



Ki(x,t) + -Q^{x,t) 



(37) 



2c s 



c s 2c[t 2 - x 2 



cl \Jclt 2 - x- 



:9(CLt-\x\) 



2cp 2 - x 2 
^/c 2 t 2 - x 2 



0(c s t-\x\) 



x 0(c s t-\x\) 
2c s t 2 ^f c lt 2 - x 2 ' 



where the kernels are: 

Ki{x,t) 

To arrive at Eq. d37l ). the prescription 9(t = 0) = 1/2 was used. The highly singular contribution 8X2 /dx in Eq. ( f37T > is 
a distribution that should be used be means of integration by parts. Proceeding as for the screw, the following dynamic PN 
equation is obtained: 



c s 0(c s t - \x\ 



(38a) 



(38b) 



dr dx' Ki(x — x', t — r)— — (x 1 , r) — — 
ox ir 



drdx'K 2 (x-x',t-r)^(x',r)-^-^(x, t) 

+ a a (x,t)=f{r 1 (x,t)). 



(39) 



Remarkably, this equation features a convolution with the sec- Whereas the term drj/dx, maximal at the dislocation center, is 
ond derivative d 2 r\/dx 2 that was not present for the screw. 
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linked to the position of the dislocation, the second-derivative 
term provides contributions from leading and trailing regions 
of the core and thus has bearings on the dislocation width. 

D. Climb edge dislocation 

Because in principle climb dislocation move diffusively, 1 
the question of their "relativistic" velocity regime might ap- 
pear irrelevant. It was however put forward some decades 
ago in a stationary context by Ang and Williams — and 
Weertman. 16 ' 42 Weertman advocated that since the atomic 
spacing across a stacking fault may differ from that across 
an unfaulted plane, an out-of-plane component of the Burg- 
ers vector, of the climb type, ought to be associated to partial 
dislocations that bound stacking faults on their glide plane. 42 
This component would then move in concert with the in-plane 
dislocation components, with the fast velocity of the latter and 
in a diffusionless fashion. Dissociation of a perfect dislocation 
into partials indeed leaves a possibility that leading and trail- 
ing partials have out-of-plane components of opposite sign, 
leading to an overall asymmetric pattern of the dissociated 
dislocation. In fact, asymmetric dissociated non-planar cores 
have recently been observed in molecular-dynamics simula- 
tions in copper, in a dynamical context where non-planarity, 
greatly enhanced by a high applied stress, arises in conjunc- 



tion with fast motion. A link to climb components might be 
that non-planar cores can be modeled by using "core field" 
corrections, 44 obtained by prescribing an effective asymmet- 
ric dipolar distribution of Burgers vectors that includes out- 
of-plane components near the core, as a model to nonlin- 
ear effects. 45 In this case however, the climb components are 
"ad hoc" and ought not be attributed the physical character 
of lattice Burgers vectors. If this connection proved true in 
the dynamical case, such corrections could be derived using 
spatial derivatives of the displacement field of climb com- 
ponents, although we shall not pursue in this direction. On 
the other hand, isolated regular climb components may en- 
ter the "dislocation part" of disconnections, a category of de- 
fects with a step associated to interfaces between grains or 
phases, which can move in a diffusionless fashion, e.g., in 
diffusionless transformation such as twinning or martensitic 
transformations. 46 Transonic twinning dislocations were re- 
ported in atomistic simulations. 47 For all these reasons, we 
find it worthwhile to write down the dynamic PN equation of 
a climb component. In the case of interfacial dislocations, and 
because we deal with isotropic media, our calculation only 
concerns situations where the phases in presence are of same 
equivalent isotropic elastic moduli (twinning, notably). 

The dynamic displacement field of a climb component is 
obtained as in the glide case with no additional complications. 
We only quote the result 



CL 



u y (*,t) = -e(t) j-j-g 



0(cLt- |x| 2 )- 



c s 2c? t 2 - |x| 



CL 



0(c s *-|x| 2 ) 



9(cLt-\x.\ 2 )- 



2c 2 t 2 



£s 

CL 



0(cst-|x| 2 ) 



9(c h t - 



X2 



xy 9(c L t- x 2 ) , Cl . , ... ... , h 

c s c^ 2 - x 2 ) y/tff - x \l c s 



(40a) 



(40b) 



Writing a yy , the resolved stress of the climb component, as: 



with the kernels 



CTyy — f-l 



„2 ( U V,V 



its value on the plane y = reads: 



^ ) ^x : x 



Vyy(x,y = 0,t) 



2 4 



(41) 



(42) 



c s 2c^ 2 



2cjt 2 - x 2 



9(c s t-\x\) 



CL yjcft 2 - X 2 

c-lx 9{c\t 



9(c L t-\x\) 
x\) 



2c s 2 * 2 JW- 
~ 9{c L t-\x\) 



(43a) 



(43b) 



The corresponding dynamic PN equation is of the form d39l ), 
where now r/(x) = u y (x,0 + ) — u y (x,0~), and where the 
coefficient of the third (instantaneous) term on the left-hand 
side (lhs) of Eq. d39l l, namely, /i/(2cs), should be replaced by 
/«Cl/(2cs) according to Eq. d42l . This case however remains 
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somewhat formal, for lack of available proper definitions of 
the associated pull-back force f'(rj). 



E. Static limit 

A first independent check of the above results consists in 
computing from Eqs. (|28]), d36al >. d36bb , (I40al > and (140bl > the 
following "static" displacement field: 



u(x) = lim / dru(x, t). 



(44) 



It is easily found that this integral applied to Eq. ( I281 l gives 
Eq. (fTST i back, and that Eqs. ( I23ab and ( I23bl i are retrieved 
with e g = l/(2cst) — > by applying it to Eqs. (1361 . Like- 



wise, the static fields [Eqs. ( I24ab and ( l24bM of the climb edge 
are retrieved from Eq. (l40l i. with the following scaling pa- 

c 2 /c 2 — 1 

rameter in the logarithm: e c = e 3 (c 1 / 2 cl/cs) l s ■ For 
both "edges", the logarithmic divergence at large sizes is re- 
placed by a divergence at large times, the true static regime 
being reached when e becomes of order the inverse system 
size. Remark in passing that the dynamic ratio e c /e g found 
here is different from its static value e c /e g = expl/(l — 
1v) = exp (c 2 /d? — l) (see Sec. Ill Cb . owing to differences 
in the limiting process employed, unless cl = e 1 / 2 cs. The 
next time-dependent correction in the asymptotic expansion 
at large times of the integral in Eq. ( f44l > is of order 0(l/t 2 ). 



where for screw dislocations, 

A(v) = yi-v*/49(l-\v\/cs\), (46a) 
B(v) = signal yV/ C 2-10(M/c s - 1), (46b) 

for glide edge dislocations (see also Ref.l26h. 



A 4 . 



Pi9(1-\v\/cl)-!£0(1-\v\/cs) 

Pl0(\v\/CL-1) 



P2 

ft 



(47a) 



B 4 

P3 °"v\/cs-l) 



x sign(w); (47b) 



and for climb edge dislocations, 

2 



4") = 2 (7)' 



in 



^(i-m/csJ-^i-h/cl) 
pi 



(48a) 



(3 2 e(\v\/c s -l) + ^6(\v\/c L -l) 



ft 



x sign(v). (48b) 



These coefficients are expressed in terms of the quantities 

Bi = |1 - (w/q) 2 ! 1 / 2 , with ci = c L , c 2 = c s , and c 3 = 
V^cs < Ci^ 



F. Stationary limit: Weertman's equations 

A less trivial independent check consists in computing the 
stationary limit of the obtained dynamic PN equations. In the 
stationary regime where the dislocation moves with constant 
velocity v, an ansatz r](x,t) = r](x — vt) should apply. It is 
observed that consistency with this ansatz requires the applied 
stress <7 a (x, t) to be either a constant, or a front moving with 
same velocity of the type a a (x, t) = a a (x~vt). Since a stress 
front necessarily propagates with one of the sound velocities, 
the dislocation velocity v is then either equal to this sound 
velocity — if the glide plane is aligned with the propagation 
direction of the front or greater — if the glide plane is inclined 
with respect to this direction. Thus, a stationary propagating 
front can only involve transonic or supersonic dislocations, 
and in this case a a (x — vt) prescribes the velocity. 

Under any of these two conditions, it is demonstrated in 
AppendixIClfor the screw and the glide edge (the climb edge is 
left to the reader) that Weertman's equations 16 are retrieved in 
the following form, which encompasses all regimes (subsonic, 
transonic for edge dislocations, and supersonic): 

- ^A(v)p.v. { dx' + (iB(v)T]'(x) 

TT J X — X 

+ a a (x) = f(ri(x)), (45) 



IV. CONCLUDING REMARKS 

To summarize, dynamic extensions of the Peierls-Nabarro 
equation were derived for screw and edge dislocations (of the 
glide and climb types) using the Green's function method pop- 
ularized by Mura, 2 and the Eshelby-type trick of using identity 
(|25| >. Besides the instantaneous term that shows up in these 
equations, the origin of which was traced to a missing distri- 
butional term in the displacements, an unexpected feature of 
the dynamic PN equations is a term involving a convolution 
with the second space derivative of the displacement jump in 
both edge cases. The obtained equations formally cover all ve- 
locity regimes, as indicated by their stationary limits. Leaving 
their solution to future work, we conclude with some remarks. 

Technically, this result was arrived at by using elementary 
Volterra solutions proportional to 5(t), which simplifies calcu- 
lations. Our expressions for the dynamic stresses induced by 
the continuous displacements r\ considerably differ from pre- 
vious results. Consider for instance a moving screw Volterra 
dislocation at time-varying position £(£), starting from rest at 
t = 0, represented by the function rj(x,t) — b0(£(i) — x) 
with £(t) = for t < 0. The time-dependent distortion 
Uz,y generated by such a dislocation has been computed by 
MarkenscofL^ 2 . Her result [also Eq. (1) of Ref. |25J consists 
of a sum of two integrals, the second one being extremely 
singular on the glide plane y — 0, added to a term that com- 
pensates for the static field of the dislocation at rest prior to 
motion. The singularity that develops in her second integral 
in the limit y — > greatly complicates the obtention of the 
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stress on the glide plane. On the contrary, this stress is readily 
deduced from our Eq. d34l i. One obtains for t > 0: 



27T(7 T 



1 



fib 



dr 



W(l-\v\) 7T 



cs Jo {t 
1 

+ - 

X 

where we wrote for brevity 



i-VT 



-v 2 c s 

;i-N) 



c s (t-r) 



(49) 



(50) 



and Un(:r,t) = v(x,t, 0). The last term in Eq. d49l stems 
from an explicit integration over times r < 0, and expresses 
the progressive erosion of the static field within a shell of ra- 
dius |x| < c s t. Note that Eq. (l49l ) remains extremely singu- 
lar: besides the Dirac term, the integral over r is ill-defined at 
r = t. As in the formalism of Markenscoff and co-workers, 
these singularities arise because the Volterra dislocation is of 
null width, and can be regularized by using smoother core 
functions. However, they do not arise in the same fashion. 
Although the physical and mathematical contents of both ap- 
proaches ought to be identical, comparing them explicitly 
proves difficult. For this reason, we found useful to give 
straightforward independent checks by detailing in the ap- 
pendix the steps leading to Weertman's equations. It should 
be noted that the latter derivation does not require separate 
consideration of the different sonic regimes, contrary to pre- 
vious works. 

Next, the instantaneous term — [/x/ (2cs)](drj/dt)(x, t), ab- 
sent from Eshelby's dynamic PN equation for screws, that we 
obtain in the dynamic equations, is of dissipative nature. It ac- 
counts for instantaneous losses by shear wave emission trans- 
verse to the slip plane as the dislocation advances. In oppo- 
sition, the nonlocal kernels represent the waves on the slip 
plane (of the shear type for screws, and of the shear and lon- 
gitudinal types for edges) that determine the core shape. In 
the subsonic steady state, transverse radiative losses in Weert- 
man's equations are exactly compensated by the energy that 
flows to the core. 26 This is the meaning of the compensation 
of terms proportional to v that occurs in the calculations of 
Appendix ICl 

Moreover, the participation of shear waves only to the in- 
stantaneous loss term of the screw and glide edges is a conse- 
quence of their in-plane character. As the instantaneous term 
in the dynamic PN equation for the climb makes clear, lon- 
gitudinal waves too would be emitted by an additional out- 
of-plane component of the Burgers vector, leading to further 
dissipation. In this connection, Gumbsch and Gac 17 already 
noted that an out-of -plane component would add some drag to 
the energetically favorable stationary radiation-free transonic 
regime 5 for glide edges that occurs at v = C3. This can be 
seen from Weertman's equations. Indeed, radiative losses are 
proportional to B(v) in Eqs. ( |46bt , J47bt and (g8b]l^ and 
whereas B(c^) = for the glide edge, this term is non-zero 
for the climb. For lack of proper knowledge about it, we pre- 
ferred not to conclude on the form of the pull-back force in the 
case of the PN equation for a climb component, but potentials 



for glide and climb edge components of the same dislocation 
ought be coupled in some manner. 

Furthermore, although admitting arbitrary velocities, it is 
not clear to us at present whether the dynamic PN equations 
should or not require the same additional regularization as 
Weertman's, their stationary limit, to produce solutions that 
behave correctly. As they stand Weertman's equations are in- 
deed known to be defective, for two different reasons. 26 First, 
in presence of a homogeneous applied stress, these equations 
admit no single-dislocation solution due to absence of dis- 
sipation in the subsonic regime v < cs where B(v) = 0. 
This problem finds its origin in the static PN equation^ To 
correct it, a phenomenological additional drag term must be 
prescribed to account for losses of lattice origin. 26 Evidently, 
such a term can be added as well in the form —a(drj/dt)(x, t) 
(a being some drag coefficient) to the lhs of our Eqs. d35b 
and d39l l, thereby "renormalizing" the already present loss 
term, but at the risk of excessively damping sounds waves. 27 
The second reason is that for a supersonic transition to take 
place, "relativistic" core contraction — a feature of the solu- 
tions to Weertman's equations that plausibly carries over to 
some extent to the instationary regime — must be forbidden 
below some microscopic scale to prevent energy from becom- 
ing infinite. 13 In the stationary limit a phenomenological im- 
plementation of this constraint consists in curing Weertman's 
equations by adding a smoothing gradient term, 26 which pro- 
vides a connection with the spatial-dispersion effects alluded 
to in the introduction. Recently, a simple device has been pro- 
posed allowing one to overcome both problems at the same 
time. It essentially consists in a suitable coarse-grained rein- 
terpretation of the PN equation in the static case. 27 A similar 
procedure could be applied to the dynamic equations. How- 
ever, it might occur that consideration of full dynamical be- 
havior alleviate the need for such regularizations. 

Finally, anisotropy has important consequences on the 
stress/velocity dependence of dislocation motion, owing to 
the presence of three sound waves in anisotropic media. 8 
The present theory could be extended to this case by appeal- 
ing to available elementary anisotropic dynamic solutions for 
displacements.— 
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Appendix A: STATIC DISPLACEMENTS BY THE GREEN'S 
FUNCTION METHOD 

This section examines only the calculation for the Volterra 
screw dislocation, as an illustration of how distributional parts 
emerge from otherwise standard Fourier integrals. For com- 
pleteness, like calculations for the glide and climb edges are 
provided in a separate document. 37 From (O and with /?*■ 
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given by Eq. ( I12ab . one has 

Tij = CijklPkl = (S i3 Sj 2 + S l2 5 j3 ) (Al) 

so that with the static Green's function [Eq. ([Totl. 



[GyfcfeT fe j](k,i) = - 



^3^2 + ^2^3 

1 



(1-1/) 



hk 3 k 2 



P* 2 (k,t). (A2) 



Then, specializing Eq. © to the static case by carrying out the 
time integration [Eq. d44b l in the first place, 



u z (x, y) = b 



dfci dk 2 e ^ klx+k2yS > k 2 

ki + ie kf + k 2 



(2n) 



(A3) 



To arrive at this integral, the FT of was carried out with the 
help of the (one-dimensional) FT of 6(—x), which evaluates 
to i/(ki + ie) with e — » + . In Eq. ( IA3b The integral over k 2 
is done first, so as to account for the prescription e — > in the 
remaining integral over ki . By contour integration, 



dk 2 k 2 e lk2V _ i 
~2tt k\ + k\ ~ 2 



sign(i/)e~ 



(A4) 



and the remaining integral over k\ is 'folded' on the positive 
semi-axis with a change of variables before letting e —> 0. 
This leads to the integral 



dfci e lkl 



ii 2n k\-\- ie 

fci 



dfci 



-fci |y I 



x 



sin(fcix) — 



■ cos(fcix) 



dfci 



sin(fcix) 
fci 



- 7T(5(fcl) 



(A5) 



The way the Dirac distribution arises in Eq. (IA5I ) makes clear 
that the prescription J +O ° dfci S(k\) = 1/2 holds. Moreover 
(G.R. 3.941-1), 

Jo h 

, . /"+ 00 dfci _ fc M . . . i ja?| 

= sign(x) / — — e W sin(fci) = sign(x) tan — , 
Jo h " \y\ 



(A6) 



so that 

r +co dfci e iklX 
7 27r fci + ie 



\ki\\v\ 



s 1 . , . Id' 
-x) H — sign(x) arctan — 

7T X 



(A7) 



Multiplying by the factor (i/2) sign(y) coming from (IA41 > 
eventually yields 

dfci dfc 2 e ^ x+k ^ fc 2 



(2tt) 2 h+ie kl + kl 

1 2/1 
= — arctan - + - sign(y)0(-x), (A8) 
Zir x 2 

whence expression ( fT6] i of u z . The edge cases are addressed 
by similar means^ 



Appendix B: DYNAMIC DISPLACEMENTS 
1. Screw dislocation 

The instantaneous screw is generated by the eigendistortion 
of nonzero component f3* y (x,t) = bS(y)0(—x)S(t). With 
now fc = (k\ + fc 2 ) 1 / 2 an d using Eq. (O, the displacement 
takes on the form 

r+°° d 3 fc 

u z [x,t) = -i J j^[G 3j k k T jk ](k,t)e lk - x 

= bc&0(t)lW(x,y,t), (Bl) 
where the following integral was introduced: 

lW(x,y,t)= f^dk 2S Mckt) kel{klX+k2y) 
J (2iry ki+ie 
_ .8 /-dfci e«*** rdk 2 sm(ctk) cik2tl 
dy J 2tt fci + ie J 2n k 

In this expression, the inner integral over k 2 is (G.R. 3.876-1): 



dk 2 sin (ctk) 



,ik 2 y 



2tt 



For ct > \y\, going to the limit e — > as in Eq. ( lA5b . the 
remaining integral is (G.R. 6.693-7): 



\jo(\ki\{c 2 t 2 - y 2 ) l l 2 )e(ct-\y\). (B3) 





f dfci 


-'J 


2tt 




dfci 


Jo 


n 


1 


1 


~2 




1 


1 



e —— Jo(N(c 2 ; 2 -2/ 2 ) 1/2 ) 
ei + ie V / 



sin(fcix) 
fci 



-irS(ki) 



Jo(fci( C 2 < 2 -y 2 ) 1/2 ) 



= -\ + - r — sin i ux ^ t2 y 2 y 1/2 ) 

2 n J u \ I 
•sign(x)0(|x| 2 -c 2 t 2 ) 



2 2 

1 fx 
I — arcsin — = 

V Vc 2 i 2 - y 2 



c 2 t 2 - Ixi: 



— arcsin 

7T 



^/ c 2 i 2 - y 



6(c 2 t 2 - |x| 2 ) 
-0(-x). (B4) 



Multiplying by (l/2)0(ct - |y|) according to Eq. dB3l . and 
differentiating the product with respect to y according to Eq. 
(lB2l yields 



xy 



%t-|x| 2 ) 



(c 2 i 2 - y 2 ) V^ 2- Nl 



+ 7rsign(j/)^(-x)d(ci- \y\) 



(B5) 



Equation (128b follows 
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2. Glide edge dislocation 

The only non-zero component of /3 is now /3* 2 (x, t) — 

b8(y)9(—x)6(t), and k k T k j — ^/3* 2 (fci<5j2 + fc 2 #ji). Thus 



G io k k T k:j =0{t)cl 



cs 



— sin(csfct)ffci(5i2 + fc2'5ii-2fcifcifc2 

/?*,. (B6) 



H su\(c\ j kt)kik\k2 

cl 



Setting fc = (k\ + fcf) 1 / 2 , the non-zero components of u are 
obtained as: 



u x (x,t) = -i I [GijkkTjk}{k,t)e 



ik-x 



, 2/)/,n f dfcidfc 2 



(2tt) 3 1 

,i(feiK+i!2l/) 



sin(csfci)fc2 
csfc(fci + ie) 



sin(cLfct) sin(csfci) 



cl 



cs 



fcl&2 
"fc 3 " 



(B7a) 



and 



u y (x, £) = — i 



+ 00 



d 3 fc 

(2^)3 



bcj6(t) 



dkidk 2 i(klX+k2v) [ sin(e s H) 
(2^) 2 1 c s fc 



sin(cLfci) sin(csfci) 



cl 



cs 



ft 2 

A- 3 



(B7b) 



In these expressions, the limit e40 was taken wherever pos- 
sible (cancellation of fci between numerator and denominator 
of fractions). Four different types of integrals are involved. 
The fist one, JW, was defined in Eq. (1B21 > and computed in 
Eq. dB31) . The three others ones are (G.R. 8.411-5 and 6.671- 
7) 



I {2) (x,y,t) = [ dkl dk2 Sin(cH) e <(*i*+faiO 
(27r) 2 fc 



(B8a) 



/ — sin(cfct)J (fc|x| : 
Jo ^ 



1 0(ct- |x| 2 ) 



2tt ^/ C 2£2 _ | x |2' 

[x,y,t)-J (27r)2 fc3 e _ flaj> 

(B8b) 

rWr-r „ rt- / dfci dfc 2 sin(cfcf)fcj <ffaa , +fari _ dJ 



(B8c) 



where the following integral was introduced: 

J(x,y,t) = -i [ dfc i dfc 2 sin(cfct)fc 2 gjfex+fc^) 

y (27r) 2 fc 3 



sign(w) 



dfci cos(fcia;) 



d^gsin((jfci|y|) 



(l + g 2 ) 3 / 2 



sin(ctfci(l + g 2 ) 1/2 ) . (B9) 



Its expression in polar coordinates shows that J is finite. The 
last equality in Eq. (IB91 > follows from elementary symmetry 
considerations and from a change of variables fc2 — > g = 
fc2/|fci|. Consider first the inner integral, and introduce for 
convenience (a and are arbitrary positive constants), 



dq qsm(bq) 
T(l + g 2 ) 3 / 2 



sin ( a(l + g 



2x1/2 



(BIO) 



This integral is not tabulated for all positive (a, b) pairs (see 
G.R. 3.875-3 for a < b). However, one integration by parts 
over q and the use of (G.R. 3.876-1) show that (Jo is the Bessel 
function) 



j{a,b) = a^-(a,b) 



dq cos(bq) 
/0 V(l + g 2 )V2 

dj , ., b 



(«(1 + ? 2 ) 1/2 ) 



= a£(a, b) + -Jo ((a 2 - 6 2 ) 1 / 2 ) 6(a - b). (Bll) 

Thus, j is a continuous solution of a homogeneous (resp. non- 
homogeneous) differential equation for a < b (resp. a > b). 
This differential equation is solved by variation of constants 
with condition j(oo, b) = (see G.R. 6.554-4 for the integra- 
tion constant). A change of variables then gives 



3{a,b) = - 



e- b -bO(a-b) 



(a 2 -6 2 ) 1/2 , 



uJq(u) du 
(u 2 + 6 2 ) 3 / 2 



It follows that 

J(x,y,t) = sign(u) 



dfci cos(fcix) 

-j(ctfci, |y|fci) 



27 sign( ^ 



lv|0(ct-|vl) 



fci 

dfci cos(fcix)e" 



-kilwl 



(c 2 * 2 -?, 2 ) 1 / 2 



i du 



(u 2 + y 2 )3/2 



x / dfci cos(fcix) Jo(ufci) 

'0 



that is, with (G.R. 3.893-2) and (G.R. 6.671-8), 
1 



(B12) 



J{x,y,t) 
- 6{ct - 



cty 
2ir 



*2 



2-2 _|_ y2 

(c 2 * 2 -*, 2 ) 1 / 2 



it du 



( w 2 + y 2)3/2 ( u 2_ 3-2)1/2 



ct y 

27T |X| 2 



l~f tX /cHi-\ x \t9(ct-\ x \ 2 ) 



(B13) 
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Integrals 1^ and 1^ follow from differentiation according 
to Eqs. dB8bt . dB8c] > as 



cf 1 



2c 2 f 2 - |x|| 



lW(x,y,t)-- 



'2tt |x|| 



cf^f 2 - |x 

x 2 |x| 2 ~c 2 t 2 (x 2 -y 2 ) 
rf^/c 2 t 2 - |x| 2 



(ct-lxk) - 2 



0(ci- |x| 2 ) 



(B14) 



Gathering all contributions within Eqs. ( |B7aK dB7bt then 
yields displacements d36a| ) and ( 136b| >. 



Appendix C: STATIONARY LIMIT 

1. Screw dislocation 

Using cr a (x, t) = <j a {x— vt) and the ansatz r)(x, t) — i](x — 
vt) (see Sec. IIII Ft in Eq. ( 127) , one sees that 77(2;) obeys the 
PN-like equation 



where 



dx' K v (x-x')r,'(x') + ^-ri'(x)+a a (x) = /'(rj(a:)), 

(CI) 



f°° P dk 

J dtK(x + vt,t) =J — e ikx K v (k) (C2) 



which features the space Fourier transform of K v (x) in the 
form of a one-sided integral over time, 



K v (k) = / dte Lkvt K{k,t). 
Jo 



(C3) 



In this expression K(k, t) is the space FT of K(x, t) [given in 
(|30l>l, which reads (G.R. 3.752-2; Ji is the Bessel function): 



K(k,t) 



ik f Cst I 

2 9(t) I Ax J c\t 2 — x 2 cos(kx) 

est Jo v 
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6{t) Ji(c s fct). 



(C4) 



The expression of K v {k) is evaluated from (IC3b with the help 
of the integrals (G.R. 6.693-1 and 6.693-2) 



dt 

T 



cos(kvt)Ji(kcst) 



= sig: 

T 



n(k)^/l-v 2 /c 2 9(l-\v\/c s ), 
sm(kvt)Ji (fccsf) 



(C5a) 



(v/cs) - sign(v)yJv 2 /4-ie(\v\/cs - 1), (C5b) 



from which: 



K v {k) = -nrsign(k)^l~v 2 /c 2 9(l - \v\/c s ) 
{v/c s ) - sign(v)^v 2 /c 2 - W(\v\/c s - 1) 



(C6) 



Since — iir sign(fc) is the FT of p. v. 1/x, the Fourier inversion 
of K v {k) is immediate as 



K v (x) = 9(1 - \v/cs\)^\-v 2 /c 2 V .v. i 
(v/c s ) - sign(u) \Jv 2 /cl - 19{\v\/c s - 1) 



S(x).(Cl) 



Putting this expression into (ICU one sees that the instanta- 
neous terms (proportional to v) cancel out mutually. Weert- 
man's equation d45b with coefficients ( 146 at and ( I46bl ) follows. 



2. Glide edge dislocation 

Again using a a (x, t) = a a (x — vt) and the ansatz rj(x, t) = 
r){x — vt) in the dynamic PN equation ( |39l ) for the glide edge, 
the resulting stationary equation takes on the form dC It where 
now 



K v (x) 



dt 



Ki(x + vt,t) 



dK 2 
dx 



(x + vt, t) 



(C8) 

in which the kernels K\ and K% are given by ( f38l >. Proceeding 
as for the screw in Sec. IC II one evaluates first the Fourier 
transforms of K\ and dK^jdx wrt. x. By means of changes 
of variable x — >• u = x/(c\J) and x —> u = x/(cst), and 
using the fact that Kx(x,t) is odd in x and that K<z(x, t) is 
even, one gets with (G.R. 3.753-5) and (G.R. 3.753-2) 

1 it 4-7 

^ l(M) = ~Yt Jl{kcst) - t sign(fc) (C9a) 



c / du 



dK- 



sin(|fc|cLttt) sin(|fc|csfu) 



2-u 2 



u 3 VT 



dx 



(k,t) = -i-kc s J (kc s t). 



— 2 ' 
(C9b) 



The next step consists in obtaining J Q °° dte lkvt Ki(k,t) in 
which we write 



ikvt 



cos(kvt) + isign(fc) sign(w) sin(|fc||«|t) (CIO) 



With (|C9at and (IClOl l. the latter integral involves the follow- 
ing integrals over time, where c stands either for Cl or for cs 
(G.R. 3.741-1,2): 

y cos(\k\vt) sxa(\k\ctu) = ^9(u- \v\/c), (Clla) 



df 



sin(|fc||u|f) sin(|/c|cfu) 



2 

4 *\u-\v\/c 



(Cllb) 



The remaining integrals over u combined with (IC lib are eval- 
uated using the pair of integrals, 
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6u~6(u-\v\/c) 2 " 



du - log 

4 



u 3 Vl — u 2 

2 „ 2 



(CI 2a) 



2-u z 



u + \v\/c 
u—\v\/cj u 3 \/l - u 2 \v\ 



f 1 du 




1, + 


/ du 




J 



log 



u+\v\/c\ 2 2-u 2 2c 1 
u—\v\/c) 4u 3 Vl - u 2 \v\u 2 



2c 2c 
|w|e \v\ 

cH 2 7T C 2 

c"6+ |v| x 2 v 2 



vT 



4u 2 



■log 



u + \v\/c 
u — \v\/c 



1 du 



1 



(«/ C )S 



= lim -j— I — - ~Jv 2 lc 2 - iO(\v\/c- 1). 



(CI 2b) 



In ( IC12bb the following transformations were applied. The 
integral is divergent at u = 0. Its divergent part is extracted 
first, and expressed in terms of x, recalling that u was intro- 
duced via the change of variables u = x/ (ct). In this form, it 
is proportional to c 2 and cancels out when assembling contri- 
butions involving cl and cs in the final step of the calculation. 



Meanwhile, the remaining finite part is integrated by parts, 
and the spurious singularity at u = \v\/c introduced by this 
transformation for \v\ < c is removed by the principal value 
prescription. Appealing next to dC5al ). ( IC5bt to deal with the 
Bessel function in (lC9a| i, these contributions to jC9a| > lead to: 



foo 2 

/ dt e lkvt K 1 (k, t) = -2in sign(fc) 4 
Jo v 



1 - v 2 /c 2 9(1 - M/c L ) + ["-!) Jl- v 2 /c 2 0(1 - H/c s 



44 



-27T sign(w)-| 



v 2 /c 2 - 1 9(\v\/c L - 1) + ( _ ^ 1 ) ^ v 2 /c 2 - 1 9(\v\/c s 1) 



Turning now to K 2 one finds, using dC9bl . dCTOt and (G.R. 6.671-7,8) 

~dK 2 



dte 



i k 1 f t 



dx 



.j , 7T . .,.(9(1 - H/c s ) . , N 7r6»(|u|/c s - 1) 
(fc, t) = -«- sign(fc) \ 1 ' 7 y + sigm»- -iL ^J_^ . 

^ v 1 - w 2 /Cs ^ v / c s ~ 1 



r 



7TU 

35- (CB) 



(C14) 



The contributions of K\, dK 2 jdx are brought back into with A(v) and B(v] given by ( 147 at and d47bb . This brings the 
K v (k), whose Fourier inversion is immediate as in the screw present edge version of ( IC1I ) down to Weertman's equation, 
case, see Eqs. iC6\ and ( IC7b . The result reads 

KJx) = A(v) p.v. - - nB(v)5(x) + ^r-S(x), (C15) 
x 2 c s 
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